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Abstract 

The classical flexure problem of non-linear incompressible elastic- 
ity is revisited assuming that the bending angle suffered by the block is 
specified instead of the usual applied moment. The general moment- 
bending angle relationship is then obtained and is shown to be de- 
pendent on only one non-dimensional parameter: the product of the 
aspect ratio of the block and the bending angle. A Maclaurin series 
expansion in this parameter is then found. The first-order term is pro- 
portional to fjL, the shear modulus of linear elasticity; the second-order 
term is identically zero, because the moment is an odd function of the 
angle; and the third-order term is proportional to /u(4/3 — 1), where j3 is 
the non- linear shear coefficient, involving third-order and fourth-order 
elasticity constants. It follows that bending experiments provide an 
alternative way of estimating this coefficient, and the results of one 
such experiment are presented. In passing, the coefficients of Rivlin's 
expansion in exact non-linear elasticity are connected to those of Lan- 
dau in weakly (fourth-order) non-linear elasticity. 

Key words: Incompressible elasticity; flexure; non-linear shear coefficient; 
experimental data. 
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1 Introduction 



The problem of flexure is one of the now classical problems of the theory of 
non-linear, incompressible elasticity. First formulated and solved by Rivlin 
PP, it has since been studied extensively in the literature, see for example 
Green and Zerna [2] and Ogden [3J. There is continuing theoretical interest 
in this problem, as can be seen, for example, in the recent study of Kanner 
and Horgan [4j. The physical problem considered is easily visualized: a 
rectangular specimen is bent by equal and opposite terminal couples applied 
on the end faces of a rectangular block, while its other faces remain free of 
traction. Rivlin pQ shows that if a circular, annular sector is assumed for 
the deformed configuration, then an elegant solution to the corresponding 
boundary value problem can be found. We recall this derivation in Sections 
ElandEJ 

Typically, the usual formulation of this boundary value problem implicitly 
assumes that the terminal moments are specified. It is shown here that 
specifying instead the bending angle through which the block is bent results 
in a simpler mathematical formulation and solution of the problem. We 
show that this solution depends on only one non-dimensional parameter: the 
product of the aspect ratio of the block and the bending angle, which we 
denote by e. If e is assumed small, then the solution at low orders has a 
particularly simple form. Because e is the product of the aspect ratio of the 
block and the bending angle, the lower-order solutions are applicable in at 
least two distinct physical regimes: the first corresponds to the bending of 
bars through an infinitesimal bending angle and the second corresponds to 
the non-linear bending of thin sheets. Although these two sub-cases are the 
most important, there are other possibilities: moderately thick blocks could 
also be considered turned through moderate angles, provided the product e 
is small. 

To analyze the problem, plane strain conditions can be assumed. Hence 
the technical interpretation of our results is that they describe the bending 
of planar sections of infinitely wide blocks. We assume that the plane strain 
results obtained are also applicable to strips with a finite, out of the plane di- 
mension larger than the block thickness, i.e. that the edge/anticlastic effects 
are negligible. This assumption is motivated by observations and measure- 
ments for the bending of rubber blocks, see Gent and Cho [5] . Of course, this 
assumption is a limitation of Rivlin's solution, because secondary fields are 
bound to be observed outside of a central area in the bent block, see Figured! 
Nonetheless, it must be kept in mind that Rivlin's solution is one of the very 
few universal solutions, valid in principle for every incompressible isotropic 
material, whatever the actual dimensions of the block and the amount of 
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Figure 1: Bending of a block of silicone rubber, with length L = 24 cm 
and thickness 2A = 2 cm. The dimension of each square drawn on the top 
surface is 1 cm x 1 cm. The picture on the right shows that even for such 
a "home-made" bending experiment, there exists a region about 3 cm wide 
around the median where plane strain is respected. 

bending. For more refined deformations, albeit limited to certain types of 
geometries, the reader is referred to the advances obtained by Shield [6]. The 
other limitation of this solution is that it might bifurcate, see [Zl El El HH] for 
an in-depth treatment of this possibility. 

Given that here, the deformation depends only on the bending angle, 
the corresponding stress distribution can be easily determined. The most 
important functional of this stress distribution is the moment that needs to 
be applied at the ends of the block to effect the deformation. In Sections H] 
and we show that the general relation between applied moment and the 
parameter e, the main relation of interest in this problem, can be expressed in 
a succinct and elegant form. It is shown that this relation is an odd function 
in e, which agrees with an intuitive expectation that the moment should be 
an odd function of the angle since the moment required to bend a block by 
an angle a say, is the opposite of the moment required to bend it by an angle 
—a. 

On expanding the moment in a Maclaurin series in e (Section EJ), we find 
that the first-order coefficient is proportional to fi, the shear modulus of 
linear elasticity. The second-order coefficient is identically zero for all elastic 
materials, because the moment is an odd function in e. This suggests that 
the linearized moment-angle relation is likely to be valid for e values beyond 
the infinitesimal range, and this is verified numerically and experimentally in 
Section |HJ Prior to this, we show in Section [7] that the third-order term in the 
expansion is proportional to /x(4 ( 5 — 1), where (3 is the non-linear coefficient of 
plane shear waves pH QH Q31 Q3] . Explicitly, 2/3 = (ji + A/2 + £>)///, where 



3 



A and T> are Landau third- and fourth-order elastic constants [15]. This 
coefficient has been measured for agar-based gels, based on the measurement 
of shear wave speeds in transient elastography p3], or on the measurement 
of homogeneous plane strain deformations [TB] . Clearly, the information 
collected from the bending of a block, such as that provided by a bending 
stiffness tester [T7], yields a simple and useful alternative to these protocols. 

In Section [5J our main results are then compared to experimental data 
obtained by performing a bending experiment on a polyurethane elastomer. 
We find that ~ 1.0. 

2 Large plane strain bending 

The fundamental assumption introduced by Rivlin [TJ to model the non-linear 
flexure of an incompressible block is that a block of length L and thickness 2 A 
is deformed under applied terminal moments into a circular, annular sector. 
For definiteness, assume that that the faces X = ±A are deformed into the 
inner and outer radii, denoted by r a , r b respectively, of the annular sector 
and the faces Y = ±L/2 are deformed into the faces 9 = ±ot, where a, the 
bending angle, is a specified constant. Plane strain conditions are assumed 
throughout. The bending angle, a, is restricted to lie in the range 

< a < 7i, (2.1) 

which only allows a block to be bent into at most a circular annulus, see 
Figure EJ 

Adopting the semi- inverse approach of Rivlin [1] , assume that 

r = r(X), 9 = 9{Y), z = Z, (2.2) 

where (X, Y, Z) and (r, 9, z) denote the Cartesian and cylindrical polar co- 
ordinates of a typical particle before and after deformation, respectively. 
Incompressibility then yields 

r 2 = 2BX + D, 9 = Y/B + C, z = Z, (2.3) 

where B, C , D, are constants. 

As noted by Rivlin [I], symmetric boundary value problems can be con- 
sidered without loss of generality and therefore C = 0. The key element in 
our solution of the bending problem is that the bending angle, a, is speci- 
fied and not the applied moment, as is usually assumed in most treatments, 
even if this assumption is implicit. It therefore follows easily that B can be 
determined as 

B = L/a. (2.4) 
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Figure 2: Sketch of Rivlin's deformation for the bending of a block made 
of an incompressible isotropic solid, with length L and thickness 2A, into a 
circular annular sector with inner and outer radii r a and r^, respectively. The 
angle a is the bending angle. The deformation is plane strain, but makes no 
assumption about the dimensions of the block or the amount of bending. 

This yields the non-homogeneous deformation field 

r = ^2(L/a)X + D, 9 = (a/L)Y, z = Z, (2.5) 

where D remains to be determined. Therefore, the inner and outer radii of 
the deformed curved surfaces are determined by 

r a , b = y/DT2(L/a)A. (2.6) 

Adding and subtracting these equations then yields 

D = (r 2 a +r 2 b )/2, r\ - r\ = AAL/a. (2.7) 

Hereafter we consider the boundary value problem where equal and op- 
posite moments are applied to the ends of the block at Y — ±L/2. The other 
classical boundary value problem of flexure, where one end is held fixed and 
a moment applied to the free end, is a subregion of the problem considered 
here with the fixed end described by 9 = 0. 

3 The stress-free boundary conditions 

The corresponding deformation gradient tensor, F, is given by 

F = diag (Ai, A 2 , A3) = diag (A, A -1 , l) , where A = L/(ar), (3.1) 
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denoting the principal stretches by Ai, A2, A3. It is now clear that Rivlin's 
solution (I2.5p is a plane strain deformation, because A3 = 1 at all times. 

For homogeneous, incompressible, elastic materials, the corresponding 
principal Cauchy stresses are given by 

T rr = -p + AiW.x, T e6 = -p + A 2 W )2 , (3.2) 

where p is an arbitrary scalar field, W = W(\i, A 2 , A3) is the strain-energy 
function and the comma subscript denotes partial differentiation with respect 
to the appropriate principal stretch. The equations of equilibrium determine 
p as 

p = J (\ x W,i - X 2 W >2 ) r _1 dr + AiW.i + K, (3.3) 
where K is an arbitrary constant. It therefore follows immediately that 

T rr = J (AiW,! - A 2 ^, 2 ) r"Mr + K, T ee = T rr + X 2 W :2 - X^. (3.4) 

Now define the function W(X) as 

W(\) = W(\,\- 1 ,1), (3.5) 

which is assumed to be a convex function. Then 

XW' = AiW,i - X 2 W j2 , (3.6) 

where the prime denotes differentiation. The stress distribution can then be 
written simply as functions of A as 

T rr = W + K, T ee = W- XW' + K. (3.7) 

The curved surfaces of the bent block are assumed to be free of traction. 
This assumption then yields 

K = -W(X a ), W(X b ) = W(X a ), (3.8) 

where 

X a = L/(ar a ), X b = L/(ar b ). (3.9) 

No assumptions have been made thus far about material symmetry. Only 
isotropic materials will be considered here. For these materials, W(Xi, A 2 , 1) = 
W(X 2l Ai, 1) , and so <^M> yields 

\ar a L J \ L ar a 

= w(±,^,i)=w(^,±M. (3.10) 
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There are two obvious solutions to these equations: r a = r&, which is physi- 
cally unacceptable, and pQ 

a 2 r a r b = L 2 , (3.11) 

which is assumed henceforth. Solving for r& and substitution into fl2.7jl o yields 
a quadratic equation for r 2 , with the following unique physically acceptable 
solution: 

^(V^+I- 2 - 4 )- (3 - i2) 

which completely determines the deformed configuration. Unusually, it is 
independent of the form of the strain-energy function (provided that (13.111) . 
which is sufficient for (I3.10P to be satisfied, is also necessary). Also note that 
for isotropic materials, it follows immediately from (I3.6P that 

W'(l) = 0. (3.13) 

Substitution of ( 13. lip and (I3.12p into (13. 9 p yields the following form for 
the stretches in terms of the non-dimensional parameter e: 



\ a = y/VT+^ + e, X b = \JVTT^-e = X~\ (3.14) 

where e is the product of the block aspect ratio by the bending angle, 

2A 

e=—a. (3.15) 
lj 

Finally, the qualitative features of the stress distribution will be deter- 
mined. Equation (13. 7j) i yields 

= ~AVP(A), = £ (2W'W + AHr'(A)) . (3.16) 

It follows immediately from these, the assumed convexity of W, and (I3.13P 
that the radial stress has a unique minimum value of —W(X a ) at r = L/a 
and since the curved surfaces are assumed stress-free, the radial stress is 
therefore compressive in the interior of the bent block. 

Convexity also yields that the hoop stress is a monotonically increasing 
function of r, compressive on the inner curved surface and tensile on the 
outer. 
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4 Some approximations of the deformed con- 
figuration 

Since the deformed configuration is independent of the form of the strain- 
energy function, asymptotic expansions in e are valid for all incompressible 
elastic materials. Expanding ( 13.1 2p and its counterpart for r\ in a Maclaurin 
series in e yields 

r a = - [1 - §e + |e 2 + G(e 3 )] , r b = ± [l + \e + §e 2 + (9(e 3 )] . (4.1) 
a a 

To the leading order, we have 

r a = r b = L/a, (4.2) 

giving a curvature 

K = a/L. (4.3) 

Thus infinitesimally thin sheets, strips, and wires are bent into a circular arc 
whose radius is inversely proportional to the bending angle. 

Truncating the series after the linear terms in (14. ip , and thus we are now 
considering thin, but not infinitesimally thin, sheets, strips, and wires, yields 

r a = L/a -A, r b = L/a + A, (4.4) 

which are again remarkably simple and which tell us that the thickness of 
the deformed block is still 2A. 

Retaining the quadratic terms in the expansions (14. ip yields that, again, 
r b -r a = 2A. 

5 Exact results for the moment 

The associated stress distribution naturally requires specification of the form 
of the strain-energy function in order to be determined. The most important 
functional of this stress distribution is the moment, M, required to bend the 
block by an angle a. This moment is given by 

M= I "Teerdr, (5.1) 

which can be re-written in terms of the parameter e defined in (I3.15P as 
follows: 

m r Xb ~ ~ 

T7T2=t- 2 W(X)X- 3 d\ + e~ 1 W(X a ), (5.2) 
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where A a and A& have been defined in terms of e through (I3.14p . To derive 
this expression, we made use of fl2.7j) o. (13. lp ?. ( 13. 7p . ()3.8P i . and (I3.15p . Note 
that in J3j p. 293], the last term of this expression is missing; note also that 
the expression of M in terms of an integral in r can be found in Rivlin's 
original paper pQ, see also Kanner and Horgan [3]. It follows trivially from 
(I3.8P 9 and ( I3.14p that M is an odd function of e. 

We remark that in this paper, M is the applied moment per unit width of 
the block: if the block's width is H, then the total applied moment is HM. 

The integral term in (I5.2p surprisingly means that explicit relations be- 
tween the moment and bending angle in terms of elementary functions are 
difficult to obtain in general. Progress can be made however for some forms of 
W. For example, Kanner and Horgan [I] compute M for the Mooney-Rivlin 
material and for the Gent [IB] material. 

We focus on the following Rivlin expansion of the strain-energy density 
in the principal invariants of the Cauchy-Green strain tensors |19j . 

W = C 01 (I-3) + C 1Q (II-3) 

+ C 02 (I - 3) 2 + C U (I - 3)(// - 3) + C 20 (II - 3) 2 , (5.3) 

where the CV,- (i+j = 1,2) are constants to be determined from experiments 
(see Erkamp et al. [16] for example), and 

/ = A 2 + A 2 + A 2 , // = A 2 A 2 + A 2 A 2 + A 2 A 2 . (5.4) 

This strain energy density is readily implemented in most finite element 
analysis packages. It includes the neo-Hookean and Mooney-Rivlin materials. 
In the case of plane strain, such as the bending problem considered here, 
I = II = X 2 + \~ 2 + 1, so that W reduces to 

W = \» [(A 2 + A^ 2 - 2) + (3(X 2 + A' 2 - 2) 2 ] , (5.5) 

where \i = 2(C i + C w ) > is the infinitesimal shear modulus, and (3 = 
2(C 2 + Cu + C 20 )//i > is the non-linear shear coefficient [HI [T2J IT3"] . 
Substitution into the general moment relation ( 15. 2 p then yields the exact, 
nonlinear relation 

JL = (1 _ 4/3 ) J e ~2 ln ( V /TT^ _ e ) + e -i v /iq^] + *p e , (5. 6 ) 

The last term in this expression shows that M grows unbounded with e, 
except in the special Mooney-Rivlin case for which (3 = and where the 
asymptotic value is M = 2[iA 2 [I]. 

In Figure El we plot M/(2fxA 2 ) as a function of e for (3 = 0.0, 0.1, 0.5, 1.0. 
A noteworthy feature of this plot is that there is a linear response quite far 
beyond the origin. This is explained in the next section. 
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Figure 3: Variations of the moment M with e, the product of the aspect 
ratio by the bending angle, for the large bending of an elastic block modeled 
by Rivlin's strain energy density ( 15. 3p . Exact results, in the cases where the 
non-linear shear coefficient (3 is equal to 0.0 (Mooney-Rivlin material), 0.1, 
0.5, and 1.0. See Figure H] for a comparison with experimental results when 
P = 1.0. 

6 Approximate results for the moment 

Expanding the general expression (15. 2p for the moment M in a Maclaurin 
series in e, and noting both that M is odd in e and (I3.13p . yields 



M 



W"(l) 



e + 



120 



W""(l) + 8W"'(1) - 3W"(1) 



e A + ... 



(6.1; 



Recall that only isotropic, incompressible materials are being considered 
here. Differentiating (13. 6p with respect to A and evaluation in the reference 
configuration yields 

W"(l) = W 1 + W 2 + W u + W 22 - 2W 12 , (6.2) 

where the partial derivatives of W are evaluated at (Ai,A2,As) = (1,1,1). 
The last of equations (6.1.88) of Ogden [3] then yields 

W"(l) = 4/x. (6.3) 

The dependence of W on A can be expressed in the form 

l^(A) = /(A 2 + A- 2 -2), (6.4) 



10 



for some function /, say. It then follows that W"(l) = 8/'(0) and W"'(l) = 
-24/'(0) = -3W"(1), or 

= -12/i. (6.5) 

So for isotropic materials the moment relation ( 16. ip can be written in the 
form 

M 



= 4^+1 W -(i)_i08^ e 3 + (9( e 



(6.6) 



^2 3 1 120 

which explains why the linear regime carries for moderate values of e in 
Figure [3J 



7 Weak non- linear elasticity 

Because we are looking at small, but not infinitesimal, elastic effects, we 
place ourselves in the theory of weak non-linear elasticity [20] ■ There, the 
strain energy density is expanded in terms of 

h = tr(E), /2 = tr(£ 2 ), J 3 = tr(£ 3 ), (7.1) 

where E is the Green-Lagrange strain tensor (with eigenvalues (A 2 — l)/2). 
For incompressible solids, Ogden [2T] shows that the expansion of W up to 
terms which are of order four or less in the Green-Lagrange strain involves 
only three material constants. In the notation of Hamilton et al. [15], it is 
written as 

W = ^I 2 + jh + Vll (7.2) 

where A and T> are nonlinear Landau elasticity constants. 

The Appendix shows that at the same order of approximation in the 
strains, the Rivlin strain energy density ( 15. 3p coincides with the fourth-order 
elasticity expansion (17. 2p when 

A = -8(C i + 2Cio), V = 2(C i + 3C 10 + 2C 02 + 2C n + 2C 20 ). (7.3) 

Conversely, we find that /3 introduced in (I5.5P can be written as 

in agreement with Zabolotzkya et al. [ITj . 

Consequently the general moment-e relation (I6.6P for fourth-order elas- 
ticity has the form 

^ 2 = l»e + l(» + A+ IVY + (9(e 5 ), (7.5) 
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or, in non-dimensional form, 



M 



§e+ |(4/3 -l)e 3 + C?(e 5 ). 



(7.6) 



fiA 2 



We remark that this expansion is in agreement with the Maclaurin series 
derived from the exact expression (15 .6p . see Rivlin [Tj. We also note that 
fourth- order elasticity is necessary and sufficient to express the onset of non- 
linearity present in the coefficient of e 3 in (16. 6ft : third-order elasticity cannot 
account for all the components of that coefficient, and we checked that fifth- 
order constants do not appear in it (these calculations are not reproduced 
here). 

Finally, we check that the expansion (17. 6p is consistent with classic elastic 
theory, which tells us that the total flexural moment required to achieve a 
curvature k of a plate with width H and thickness 2A is 



where E is Young's modulus and / is the moment of inertia of cross-sectional 
area. Here the block is rectangular and / = H(2A) 3 /12. Also, the curvature 
is k = a/L, see Section HI Moreover, it is known that in classic plane strain 
theory, E = 2fi(l + v)/(l — u 2 ), where v is Poisson's ratio. As we are dealing 
only with incompressible solids, v = 1/2, and we find that the linear term in 
(17.61) . or in (16. 6 p (or the term obtained from a linear expansion in e of (15. 6p ) 
is indeed given by (17. 7p . 

8 Experimental results 

We conducted bending experiments on several strips of elastomers, using a 
Tinius Olsen bending stiffness tester. That tester meets the requirements 
of the ASTM standard test method E855 [TTJ . We used strips which were 
about 4.5 mm thick. The tester imposes a moment at two points of the 
strip separated by half-an-inch. Hence the aspect ratio of the strips, with 
respect to the bending experiments, was A/L = 4.5/12.7 ~ 0.35. We bent 
the samples by small, moderate, and large bending angles, but noticed that 
at large angles, pinching (and perhaps also wrinkling |8J) took place on the 
inner face of the bent strip. Consequently we only retained the data up to 
an angle of 60°, which gives the following range for the expansion parameter: 
< e < 0.37. 

In Figure HJ we show the results of one representative experiment, for a 
strip of polyurethane elastomer, shore hardness 40A. On the vertical axis, the 



HM 



EIk 



(7.7) 
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Figure 4: Bending of a strip of polyurethane: variation of the non- 
dimensional measure of the bending moment M/M m with e, the product 
of the strip's aspect ratio by the angle of bending. Circles: experimental 
data; Dashed straight line: linear elasticity, equation (17. 7ft : Full line plot: fit- 
ting of third- and fourth-order elasticity effects with the data, by adjusting 
the non-linear parameter to the value 1.0 in (17.61) . 

variable is a non-dimensional measure of the moment, M/M m , where M m is 
the maximum bending moment of the tester's pendulum (see ASTM standard 
test method E855 [17J for details); the actual value of M m is irrelevant, as 
we are only interested in measuring the non-dimensional parameter 0. On 
the horizontal axis, the variable is e. The circles represent the recorded 
experimental data (16 measurements in the 0-60° degree range for the angle 
of bending). The straight dashed line corresponds to the fitting with linear 
elasticity theory (/3 = 0, equation (17.71) ). and the full thick plot corresponds to 
the fitting with fourth-order elasticity theory (17. 6p . Only one parameter ((3) 
is to be determined from the cubic relation ( 17. 6p . which ensures the existence 
and unicity of the fitting parameter f3. We obtained a good agreement when 
/3 = 1.0. 
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Table 1. Experimental results collected for the bending of a polyurethane 
strip with aspect ratio A/L ~ 0.35. The first line gives the bending angle a 
in degrees; the second line gives a measure of the moment M , up to a 

multiplicative factor. 
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APPENDIX: Correspondence between exact 
non-linear elasticity and weakly non-linear elas- 
ticity 

In the exact (finite) theory of non-linear elasticity, there are no restrictions 
to impose on the magnitude of the strain. Often the strain-energy density W 
is written in terms of the first three principal invariants of the Cauchy-Green 
right strain tensor C, 

J = tr(C), // = |[(trC) 2 -tr(C 2 )], 7i7 = detC. (A-l) 

For incompressible solids, III = 1 at all times, and W = W(I, II) only. 

In the weakly non-linear theory of elasticity, W is expanded up to a certain 
order in a certain measure of strain, and all higher order terms are neglected. 
Often the Green-Lagrange strain tensor E = (C — I)/2 is favoured, and the 
expansion is made in terms of the quantities I\, I2, h defined in ( 17. ip . 

There exist of course connections between the two theories. For instance, 
Rivlin and Saunders [22] show that the Mooney strain-energy density of exact 
non-linear incompressible elasticity, 

W = C al {I-3) + C 10 (II-3), (A-2) 
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coincides, at the same order of approximation, with the general weakly non- 
linear third-order elasticity expansion, 

W = fil 2 + (A-3) 

The connections between the material constants Cqi, Cio, /i, and A are 

H = 2(C 01 + C W ), A = -8(C 01 + 2C W ), (A-4) 



or conversely, 



Coi = i(2/i+j), C 10 = - 1 -L + y ), (A-5) 



see Goriely et al. [23], where there is a misprint. Now we show how Rivlin's 
strain-energy density ( 15. 3p is connected to the fourth-order elasticity expan- 
sion ( 17. 2p for incompressible solids. 

The general relations between J, II, III, and Ii, I 2 , I3 are well-known 
and straight-forward to derive: 

1 = 3 + 2I U 
II = 3 + Ah- 2I 2 + 2I\, 
111 = 1 + 2h + 2I\ - 2I 2 + pf - Ahh + |J 3 . (A-6) 

For incompressible solids, III = 1 is enforced at all times, and so 

h = -II +I 2 - pf + 2hh - |/ 3 , (A-7) 

showing that I\ is at least of order 2. Squaring gives 

I\ = I\ + H.O.T., (A-8) 

where "H.O.T." is the acronym for "Higher Order Terms" (here, higher than 
fourth-order terms). Multiplying f lX^7|) by I 2 yields IJ 2 = I 2 2 -Ifl 2 + H.O.T. 
or, using (1A-8j) . 

hl 2 = II + H.O.T. (A-9) 
Substituting (TOP and ([OD into (lA"^7j) gives [151 EH 

h = h-\h + 11 + H.O.T.. (A- 10) 

Hence, the relations (1A-6|) reduce, for incompressible solids, to 

/ - 3 = 2I 2 - §J 3 + 2/f + H.O.T., 
// - 3 = 2I 2 - fl 3 + 6/| + H.O.T., 

(/ - 3) 2 = (J - 3)(/J - 3) = (// - 3) 2 = All + H.O.T., (A-ll) 
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and Rivlin's strain-energy density ( 15. 3 p reduces to 

W = 2(C 01 + C w )h - (C 01 + C 10 )h + 2(C 01 + 3C 10 + 2C 20 + 2C n + 2C 02 )/ 2 2 . 

(A-12) 

Clearly, it coincides with the fourth-order elasticity expansion (17. 2p , with the 
connections ( 17. 3p . 
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